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Near-surface effects and solar-age determination 
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The dominant part of the difference between the observed and model frequencies of the Sun can be approximated by 
a power law. We show that when this empirical law is employed to correct the model frequencies and then the small 
frequency separations, Svo2(n), are used for solar age determination, the results are consistent with the meteoritic age 
(4.563 Gyr < t < 4.576 Gyr). We present the results and compare with those obtained by using the ratios, ro2(n), of small 
to large frequency separations. 
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1 Introduction 

It is known that there is a systematic offset between the ob- 
served and model frequencies of the Sun. This offset in- 
creases with increasing frequency and is shown (Kjeldsen 
et al. 2008) to be fitted well with a power law as 



^obs(") - ^bost(^) = a 



^obs(™) 



(1) 



where the power b is determined to be 4.90 for the Sun. 
Here f hs( n ), and i^b ea t(ri), represent the observed, and 
the best model, frequencies with spherical degree 1=0, and 
radial order n; vq is a constant frequency (chosen to be 
3100 /iHz for the Sun), which corresponds to the frequency 
for peak power in the spectrum. This difference between 
the observed and calculated frequencies exists due to im- 
proper modelling of the outer turbulent convective layers of 
the Sun. The outer layers affect the high frequencies most, 
as the upper turning point of the high-frequency waves 
are closer to the surface. Since all stellar models are cal- 
ibrated with respect to the Sun, and since, thanks to the 
recent developments, we are at a stage to have observa- 
tions of individual frequencies of stars other than the Sun, 
it is important to understand the effects of near-surface 
stellar layers on the oscillation frequencies. These effects 
also influence the small frequency separations, 8vi 1+2 (n) = 
v n \ — Vn-i 1+2, and hence, for instance, inferences of stel- 
lar ages. Here v n \ is the frequency of a mode with spher- 
ical degree I and radial order n. The use of 8vn + 2(n) for 
the purpose of solar-age determination has been presented 
in earlier works (e.g. Dziembowski et al. 1999; Bonanno 
et al. 2002). It has also been shown by several authors 
(e.g. Roxburgh & Vorontsov 2003; Otf Floranes et al. 2005; 
Christensen-Dalsgaard 2009) that the frequency separation 
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ratios, r u+2 {n) = {v n i-v n ^u +2 ) / {vni-v n -u), are not as 
sensitive to the near-surface layers. Christensen-Dalsgaard 
(2009) provides a comparison between the use of 5vi i +2 (n) 
and ru + 2(n). Here, we present our preliminary results on 
comparing different ways of using the seismic data for de- 
termination of the solar age, including the application of 
near-surface correction. 



2 Tools and methods 

We computed a series of standard solar models, using 
the stellar evolution codes, ASTEC (Christensen-Dalsgaard 
2008a) and GARSTEC (Weiss & Schlattl 2008), with dif- 
ferent ages while keeping the luminosity, and the surface 
value of Z/X fixed to 3.846 x 10 33 erg s" 1 , and 0.0245 
(Grevesse & Noels 1993), respectively. Here Z is the mass 
fraction of the elements heavier than helium, and X is 
that of the hydrogen. Models were computed with the 
OPAL equation of state (Rogers & Nayfonov 2002), OPAL 
opacities (Iglesias & Rogers 1996) together with the low- 
temperature opacities from Alexander & Ferguson 1994, 
and using the Adelberger et al. (1998) or NACRE (Angulo 
et al. 1999) nuclear reaction rates. We considered both the 
commonly used value R\ = 6.9599 x 10 10 cm (Auwers 
1 891) of the solar radius, and the value R 2 = 6.9551 x 10 10 
cm found by Brown & Christensen-Dalsgaard (1998). Mod- 
els were computed both starting at the zero age main se- 
quence (ZAMS) and including pre-main-sequence (PMS) 
evolution. We calculated the frequencies of our models us- 
ing ADIPLS (Christensen-Dalsgaard 2008b), and compared 
the small frequency separation, 5vo 2 {n), and ro2(n) of the 
models with those of the Sun (BISON, Chaplin et al. 2007). 
We selected the best model minimizing the following \ 2 , 
adapted for the separation ratios, r 02 (n), when relevant: 
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Table 1 Ages of the best models 





i?i =6.9599xl0 10 cm 


R 2 = 6.9551 xl0 10 cm 


Evolution code + 


Seismic property 


Age* (before 


Age* (after 


Age* (before 


Age* (after 


Nuc. reaction rates: 


used in x 2 


correction)(Gyr) 


correction)(Gyr) 


correction)(Gyr) 


correction)(Gyr) 


ASTEC + 


5fo2 (n) 


4.58±0.05 


4.56±0.05 


4.60±0.06 


4.57±0.06 


Angulo etal. (1999) 


ro2(n) 


4.55±0.04 


4.55±0.04 


4.55±0.05 


4.55±0.05 


ASTEC + 


5v 02 {n) 


4.57±0.05 


4.55±0.05 


4.59±0.05 


4.56±0.06 


Adelberger et al. (1998) 


ro2{n) 


4.54±0.05 


4.54±0.05 


4.54±0.05 


4.54±0.05 


GARSTEC + 


5voi(n) 


4.60±0.08 


4.56±0.08 


4.61±0.08 

(PMS: 4.67±0.08) 


4.56±0.08 

(PMS: 4.62±0.08) 


Angulo etal. (1999) 


ro2(n) 


4.57±0.08 


4.57±0.08 


4.57±0.08 

(PMS: 4.62±0.08) 


4.57±0.08 

(PMS: 4.62±0.08) 



* Evolution starts with ZAMS (Zero Age Main Sequence) unless otherwise stated. 



3 Results and discussion 



X 2 (M) 2 ) 



Y 



X,, mo del 
OU 02 



(n)} 2 



o[8vQ 2 {n)] 2 



where N is the number of modes, Sv^i 71 )^ an d ^02° del ( n ) 
are the small frequency separations from the observations 
and the models; a[Si , 02(n)] represents the uncertainty in the 
observed Svo2(n). We repeated this comparison for 5v 2{n) 
and r 2(n) of the surface-corrected frequencies (cf. Eq. 1). 




4.45 



4.50 4.55 4.60 
age(Gyr) 



4.65 4.70 



Fig.l The results with R lt using NACRE (Angulo et 
al. 1999). We present the x 2 obtained using <5fo2(«) of the 
uncorrected model frequencies (solid line), Svo2 (n) of the 
corrected model frequencies (dashed line), r 2(n) of the 
uncorrected model frequencies (squares), and r 02 (n) of the 
corrected model frequencies (crosses). 



The results of the age determinations using different input, 
as described in Section 2, are summarized in Table 1, and 
one of the cases is illustrated in Fig. 1. 

The ages of our best models found using the two dif- 
ferent evolution codes are consistent within the error limits, 
differing by 0.7% at most, and they are also compatible with 
the values in the literature obtained by employing small fre- 
quency separations: 4.57 ± 0. 1 1 Gyr (Bonanno et al. 2002), 
and 4.66 ±0.11 Gyr (Dziembowski et al. 1999). Moreover, 
our two sets of results agree with the meteoritic age, 4.563 
Gyr < t < 4.576 Gyr (Wasserburg, in Bahcall & Pinson- 
neault 1995), within the uncertainty limits. 

In principle we can determine the duration of the PMS 
evolution, up to the ZAMS, by calibrating the models using 
PMS evolution to get the meteoritic age of the Sun (see, 
e.g. Morel et al. 2000, for a discussion on the definition 
of ZAMS). Our preliminary analysis suggests that the PMS 
evolution for the Sun is 0.05±0.1 1 Gyr in order to be con- 
sistent with the meteoritic age. A more precise and detailed 
analysis is obviously needed for defining ZAMS. 

The results obtained using separation ratios are consis- 
tent with those obtained using the small frequency separa- 
tion of the corrected frequencies. However, since the appli- 
cation of the surface correction has no significant effect on 
the age determined by the separation ratios, using these ra- 
tios yields a more robust age determination. 
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